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Understanding the interaction between cavity photons and electronic nanocircuits is crucial for 
the development of Mesoscopic Quantum Electrodynamics (QED). One has to combine ingredients 
from atomic Cavity QED, like orbital degrees of freedom, with tunneling physics and strong cavity 
held inhomogeneities, specihc to superconducting circuit QED. It is therefore necessary to introduce 
a formalism which bridges between these two domains. We develop a general method based on a 
photonic pseudo-potential to describe the electric coupling between electrons in a nanocircuit and 
cavity photons. In this picture, photons can induce simultaneously orbital energy shifts, tunneling, 
and local orbital transitions. We study in details the elementary example of a single quantum dot 
with a single normal metal reservoir, coupled to a cavity. Photon-induced tunneling terms lead to a 
non-universal relation between the cavity frequency pull and the damping pull. Our formalism can 
also be applied to multi quantum dot circuits, molecular circuits, quantum point contacts, metallic 
tunnel junctions, and superconducting nanostructures enclosing Andreev bound states or Majorana 
bound states, for instance. 


I. INTRODUCTION 


Cavity Quantum ElectroDynamics (Cavity QED) en¬ 
ables the study of the interaction between light and mat¬ 
ter at the most elementary level, thanks to the achieve¬ 
ment of a strong coupling between a single atom and a 
single photon trapped in a microwave or optical cavity^. 
This paradigm has been recently brought into supercon¬ 
ducting circuits: artificial atoms consisting of two level 
superconducting circuits have been coupled to super¬ 
conducting cavities^, in the context of Circuit QED. 
These experiments provide an ideal playground to test 
the basic laws of quantum mechanics because they can 
be described in terms of simple models like the Jaynes- 
Cummings Hamiltonian. However, such models conceal 
essential physical differences between Cavity and Circuit 
QED. On the one hand, the coupling between isolated 
atoms and cavity photons mainly occurs due to the sen¬ 
sitivity of the atom electric dipole to the cavity electric 
field. This coupling depends on microscopic details since 
the atomic dipole is set by the structure of the atom elec¬ 
tronic orbitals. Furthermore, one can generally perform 
the ’’electric-dipole approximation” which assumes that 
the cavity field varies little on the scale of the atomic 
system^. On the other hand, the behavior of submi- 
cronic superconducting circuits is essentially insensitive 
to microscopic details due to the rigidity of the super¬ 
conducting phased. For instance, the behavior of a su¬ 
perconducting charge qubit can be described with one 
macroscopic variable, i.e. the total charge of a super¬ 
conducting island^i. This charge can vary due to the 
Josephson coupling between the island and an external 
superconducting reservoir. The coupling between the su¬ 
perconducting charge qubit and the cavity is usually de¬ 


scribed in terms of a capacitive coupling between the 
superconducting island and the cavity central conductor. 
As a result, the chemical potential of the superconduct¬ 
ing island is shifted proportionally to the cavity electric 
potential^. This picture implies strong inhomogeneities 
of the photonic electric field on the scale of the supercon¬ 
ducting qubit, in contrast to what is generally considered 
in atomic Cavity QED for a single atom. 

Recent technological progress is enabling the develop¬ 
ment of a new type of experiments where nanocircuits 
based on carbon nanotubes, semiconducting nanowires, 
two-dimensional electron gases or graphene, are cou¬ 
pled to coplanar microwave cavities^Ti^S. This paves 
the way for the development of ’’Mesoscopic QED”, 
a denomination introduced in a pioneering theory 
worb^i. Mesoscopic QED opens many possibilities be¬ 
cause nanoconductors can be tunnel-coupled to vari¬ 
ous types of fermionic reservoirs such as normal met¬ 
als, ferromagnets^^ or superconductors^^, in a large 
variety of geometries. So far, theoretical studies on 
Mesoscopic QED have mainly focused on quantum dot 
circuits^i2i“— . Several configurations have been sug¬ 
gested to reach the strong coupling regime between an 
electronic spin and cavity photonsSir— , or more gener¬ 
ally, to develop quantum computing schemes^^— . Meso¬ 
scopic QED also tackles problems which go beyond the 
mechanics of closed two level systems coupled to cavities, 
usually considered in Cavity or Circuit QED. The inter¬ 
action between electronic transport and the light-matter 
interaction leads to a rich phenomenology2ii^“ii2ii^i^. 
Besides, coplanar cavities could be used as a power¬ 
ful probe to reveal some exotic properties of hybrid 
nanocircuits, like for instance the existence of topolog¬ 
ical superconducting phases^, Majorana quasiparticle 
modes^^— or spin-entanglement in a Cooper pair 
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beam splitter— On the experimental side, pioneer¬ 
ing works have focused on mesoscopic rings^ and metal¬ 
lic tunnel junctions^. More recently, experiments have 
been performed with single quantum dots^iiS and double 
quantum dots (DQDs)ii“— with normal metal reservoirs. 
Reaching the strong coupling regime between the charge 
states of a DQD and a cavity is still a challenge^. Nev¬ 
ertheless, interesting resonance phenomena have already 
been observedii^—di. Several experiments have also pro¬ 
vided evidence for a modification of the cavity behavior 
by hnite bias transport through a DQDiiiiS, including a 
maser effect^S. 

These recent developments call for a full description 
of the coupling between a hybrid nanocircuit and cavity 
photons. One question naturally arises: is Mesoscopic 
QED closer to atomic Cavity QED or superconducting 
Circuit QED? What are the specificities of the coupling 
between a nanocircuit and a cavity? So far, most theory 
works have considered a capacitive coupling between the 
nanocircuit and the cavity central conductor, by analogy 
with Circuit QEDSL™i^i^-i0i43^^_ approach im¬ 
plies a coarse grained electric description of the nanocir¬ 
cuit, and a concentration of the non-homogeneous pho¬ 
tonic electric field inside some capacitive elements. A few 
works have considered a direct coupling between the mo¬ 
tion of electrons trapped in the nanoconductors and the 
bare cavity electric field, which is assumed to be constant 
on the scale of the nanocircuit^^— This is remi¬ 
niscent of the descriptions used in Cavity QED^. In this 
article, we introduce a description of Mesoscopic QED 
which bridges between these two approaches. We use 
a model which focuses on conduction electrons tunneling 
between the different elements of a nanocircuit. The tun¬ 
neling electrons occupy quasi-localized orbitals in each 
nanocircuit element, which recalls the atomic orbital de¬ 
gree of freedom of Cavity QED. However, there also exists 
collective plasmonic modes in the nanocircuit, which can 
screen at least partially the cavity fields. We use a gauge- 
invariant Mesoscopic QED Hamiltonian which accounts 
for the non-uniform screening of the cavity fields inside 
the nanocircuit, and for the electromagnetic boundary 
conditions provided by cavity conductors and voltage- 
biased nanocircuit DC gates. In the limit where photon- 
induced magnetic effects are negligible, we can reexpress 
the Mesoscopic QED Hamiltonian in terms of a scalar 
photonic pseudo-potential. This picture unifies the differ¬ 
ent approaches used so far to describe Mesoscopic QED 
devices, since the photonic pseudo-potential can vary lin¬ 
early with space in the case of a locally uniform photonic 
electric field (dipolar coupling limit), as well as remain 
constant inside a given circuit element in the limit of 
a coarse grained circuit model. In the framework of a 
tunneling model, the photonic pseudo-potential leads to 
various types of linear electron/photon coupling terms: 
cavity photons shift the orbital energies of the differ¬ 
ent nanocircuit elements, but also induce simultaneously 
tunneling and local orbital transitions. This general de¬ 
scription can be used to study the behavior of many dif¬ 


ferent types of Mesoscopic QED devices. For instance, 
it can be applied to quantum dot circuits, molecular cir¬ 
cuits, quantum point contacts, metallic tunnel junctions, 
and superconducting nanostructures enclosing Andreev 
bound states or Majorana bound states. To illustrate 
the richness of our approach, we study in details the el¬ 
ementary example of a cavity coupled to a ’’quantum 
RC circuit”, i.e. a single quantum dot coupled to a sin¬ 
gle normal metal reservoir. The photon-induced tunnel¬ 
ing terms between the quantum dot and the reservoir 
induce a non-universal relation between the cavity fre¬ 
quency pull and the cavity damping pull, contrarily to 
what is expected with purely capacitive coupling schemes 
at low temperatures^^— . 

This paper is organized as follows. In section H A, 
we discuss the gauge-invariant Mesoscopic QED Hamilto¬ 
nian, which involves a photonic vector potential. In sec¬ 
tion H B, we perform a unitary transformation to obtain 
a new Hamiltonian where the electron/photon coupling 
is due to the scalar photonic pseudo-potential. In section 
H C, we reexpress the photonic pseudo-potential scheme 
in the framework of a tunnelling model. In section H D, 
we discuss the application of our formalism to the case of 
nanocircuits with superconducting elements. In particu¬ 
lar, we give an explicit Hamiltonian for a nanostructure 
with Majorana bound states coupled to a cavity. In sec¬ 
tion HI, we work out in details the case of the quantum 
RC circuit coupled to a cavity. Section IV concludes. 
Appendix A gives a detailed mathematical justification 
for the form of the Hamiltonian of section H A, on the 
basis of an effective model which separates physically the 
tunneling electrons occupying individual orbital states in 
the different elements of a nanocircuit, from the plas¬ 
monic collective modes of this nanocircuit. Appendix B 
discusses the advantages of the photonic pseudo-potential 
scheme. 


II. GENERAL DESCRIPTION OF 
MESOSCOPIC QED 

A. Gauge invariance and minimal conpling scheme 

A nanocircuit encloses a large variety of degrees of free¬ 
dom. Our main interest in this paper is the interaction of 
the cavity with the ensemble T of the ’’tunneling” con¬ 
duction electrons, which can occupy orbital states in the 
different elements of a nanocircuit and tunnel between 
them. However, there also exists plasmonic electronic 
modes which are responsive for the screening of exter¬ 
nal fields from a massive metallic electrode. Plasmons 
can also convey screening currents in the reservoir elec¬ 
trical lines, which reequilibrate locally charge distribu¬ 
tions on tunnel junctions after a tunneling event. The 
question of plasmonic modes in a nanocircuit is compli¬ 
cated since nanoconductors have generally a low elec¬ 
tronic density, which allows only a partial screening of 
the cavity field^. In principle, even the metallic con- 
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cavity conductor 



FIG. 1: Scheme of a loopless nanocircuit with source/drain 
electrodes (blue) and electrostatic gates (green) embedded in 
a photonic cavity (purple). The yellow cloud represents the 
photonic field. The cavity presents some protrusions (pur¬ 
ple stripes) fabricated to increase the photonic field inhomo¬ 
geneities (darker yellow areas). 


tacts of the nanocircuit are not able to screen totally the 
cavity field due to their thinness. In the limit where the 
nanocircuit contacts are connected to outside electrical 
lines with a low impedance, it is reasonable to assume 
that the dynamics of plasmonic modes in the nanocir¬ 
cuit is very fast compared to the dynamics of the cavity 
modes and of the tunneling electrons^. In this case, one 
can assume that the plasmonic modes in the nanocircuit 
generate a distribution of screening charges which is en¬ 
slaved to the position of the tunneling charges and to the 
value of the cavity field. This produces a renormaliza¬ 
tion of the cavity field properties, and in particular its 
spatial dependence nearby and inside the nanocircuit. It 
also produces a renormalization of the Coulomb inter¬ 
actions between the tunneling charges. In general, the 
Coulomb interactions between the nanocircuit charges is 
also renormalized by the screening charges on the cav¬ 
ity conductors (see Appendix A). The positive ions and 
valence electrons constituting the crystalline structure of 
the nanocircuit can be treated in a mean field approach, 
leading to an effective potential VconfC^) which confines 
tunneling electrons inside the nanocircuit. This confine¬ 
ment potential can delimit quantum dots, or induce dis¬ 
order effects in nanoconductors, for instance. This ap¬ 
proach is sufficient provided the background charges pro¬ 
ducing Vconfi'^) and the microwave cavity are far off res¬ 
onant. In this framework, gauge invariance imposes the 
Hamiltonian 


Htot = j (r)/ir(r)'0(rO + (1) 


with 

C -^ ' _^ \ 2 

-I- eA(r) I -e$/jarm(^)-el4on/C^) , 

( 2 ) 

Hcoul > 

( 3 ) 

A(7=) =4f)f(o-at) , (4) 

and e > 0 the elementary charge. The field operator 
includes all the tunneling charges of the nanocircuit. 
Above, h'j-{r) is a single electron Hamiltonian, G{r,r') 
describes the Coulomb interaction between the tunneling 
electrons (see Appendix A for details), and hujQd^d is the 
Hamiltonian of a (renormalized) cavity electromagnetic 
mode. The potential ^harmC^) is due to DC voltages 
applied on the nanocircuit electrostatic gates and to con¬ 
stant charges on floating conductors like the cavity cen¬ 
tral conductor in a coplanar geometry. The coupling be¬ 
tween the cavity and the tunneling electrons arises from 
the photonic vector potential term of Eq.ljT]). For sim¬ 
plicity, we have considered a single cavity mode with cre¬ 
ation operator d^ and vector potential profile A{r). This 
description can be generalized straightforwardly to the 
multimode case by introducing a sum on a cavity mode 
index in Eqs. (HD and (HD- Note that A(r) can strongly 
vary on the scale of the nanocircuit. 

The formal description of the electromagnetic field and 
the light/matter interaction in Mesoscopic QED requires 
some care. In atomic cavity QED, the cavity field is gen¬ 
erally assumed to vary slowly on the scale of an atom^. In 
contrast, we have to take into account strong spatial vari¬ 
ations of the photonic field on the scale of a nanocircuit. 

In particular, a protrusion of a cavity conductor can be 
used to increase the photonic field locally, close to a given 
quantum do t^^d^d^d^ (See Fig. I). The DC voltage- 
biased electrostatic gates, used to control the nanocircuit 
spectrum, provide supplementary boundary conditions 
on the electromagnetic field. Besides, we have already 
mentioned that the plasmonic screening charges on the 
nanocircuit conductors can modify the photonic vector 
potential profile, especially inside and around fermionic 
reservoirs. In order to describe this complex reality and 
justify mathematically the shape of Eqs. (HMD , we intro¬ 
duce in Appendix A an effective model which separates 
physically the tunneling electrons occupying individual 
orbital states in the nanocircuit from the collective plas¬ 
monic modes of the nanocircuit. In this framework, one 
can introduce a Hodge decomposition of the electromag¬ 
netic fiel d^^'^° , which can be obtained thanks to the use 
of an auxiliary el^trostatic Green’s function. This leads 

to the terms in A(r), G{r,r') and ^harm(^) in Eqs. 

(HMD- 

In this section, we disregard the spin of charge carri¬ 
ers because we want to focus on regimes where charge 
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effects prevail. The introduction of the spin degree of 
freedom in Hamiltonian ([1]) would be straightforward by 
invoking gauge invariance (see for instance section PlE 11) . 
The direct coupling between the particles spin and the 
photonic magnetic field can usually be disregarded be¬ 
cause this coupling is expected to be very weak unless 
collective excitations are considered in a large ensem¬ 
ble of spins^ir— . From recent predictions, realSl^— or 
effective22r— spin-orbit coupling provide an alternative 
way to obtain a spin-photon coupling. 

In Bp.®, the coupling of tunneling electrons to cavity 

photons occurs through two terms, one in H -|- H , 
corresponding to single photon transitions, and the other 
one in corresponding to two photon transitions. The 
second has no reason to be disregarded in the general 
case, even if it has no structure in the electronic sec¬ 
tor. For instance, this term is crucial for determining 
the existence of a Dicke quantum phase transition when 
many two-level systems are coupled transversely to a 
cavity2iiS. The term can also bring a non-negligible 
contribution to the cavity frequency pull caused by one 
nanocircuit, which is a central quantity in Mesoscopic 
QED experiments (see Appendix B). Hence, it is conve¬ 
nient to introduce a new representation in which the A^ 
term is eliminated. This task is completed in the next 
section. 


B. Photonic pseudo-potential scheme 


In atomic Cavity QED, the Power-Zienau-Wooley 
(PZW) transformation formally enables the elimina¬ 
tion of a space-dependent A(f^^ term from the sys¬ 
tem Hamiltonian^. This generalizes the widely known 
’’electric-dipole” transformation used in the case where 
A(fJ can be considered as constant on the scale of an 
atom. However, the natural variables of the PZW Hamil¬ 
tonian are polarization and magnetization densities as¬ 
sociated to the charge distribution of an atomic system, 
which are not directly useful in our case. To reexpress 
the PZW Hamiltonian in terms of particles coordinates 
and momenta, it is necessary to perform a multipolar 
development of the light/matter interaction^. In prac¬ 
tice, this development is performed to a limited order, 
which means that only a moderate space dependence of 
the photonic electric and magnetic fields is taken into 
account. In this section, we perform a unitary trans¬ 
formation of Hamiltonian © in the same spirit as the 
PZW transformation, but with specificities required for 
Mesoscopic QED. Upon disregarding magnetic effects, we 
can take into account a strong space dependence of the 
photonic electric field, thanks to the use of a scalar pho¬ 
tonic pseudo-potential. This procedure avoids the de¬ 
velopment of the light/matter interaction in terms of 
multipoles. This is necessary to describe e.g. multi¬ 
quantum dot circuits with strongly asymmetric capac¬ 
itive couplings to the cavity. Furthermore, we use a 


second-quantized description of charges, which is more 
convenient to describe central phenomena of mesoscopic 
physics like tunneling but also superconductivity (see 
next section). 

To account for the possibly strong spatial variations 
of A{r} on the scale of a nanocircuit, we define inside 
the volume of the nanocircuit a reference point tq and a 
continuous functional path C{f,f') relating any points 
r and r'. We also define a photonic pseudo-potential 
V±{l^){a + a'^) with 

V^{^)=u;o[ 1(1^').fe. (5) 

J C(fo,f) 

Importantly, one can use Vf.Vj_(l^) ~ u:QA{r) inside the 
nanocircuit volume, provided the flux of Vf? A .4 through 
the nanocircuit can be disregarded. This criterion is 
generally satisfied, considering the typical amplitude of 
the cavity magnetic field and the size of nanoconductors. 
In this case, the arbitrariness on the choice of the con¬ 
tours C{r, r ') leads only to marginal effects. One excep¬ 
tion is nanocircuits with large loops leading to magnetic 
Aharonov-Bohm effects. This case is beyond the scope of 
the rest of this article, which focuses on photon-induced 
electric effects. We perform the transformation 

SZt = (6) 


with 


U = exp 


e(a — d’f) 
hioo 


J (fir . (7) 


This gives 

(firtp'^ {r)hr{r)tl>{r) + Hcoui + fu^oo^a 
+ V{a +a)) + /hwA (8) 

with 

V = —e J (firV±{'fi')tp^ {r)tp{r} (9) 

and hr{r) = -h?/2m- - eVconffifi)- In 

this new representation, the single electron potential and 
the Coulomb interactions are renormalized by the term 
in V^, but the electron/photon coupling takes a simpler 
linear form. The motion of tunneling electrons is modi¬ 
fied by the pseudo-potential V±{'fi){a + d’^). Note that 
the arbitrariness on the choice of ro is not an issue, since 
a change in rg shifts VI ("t^) by a global constant and 

leads to a unitary transformation of Htot- In tbe limit 
where the spatial dependence of the photonic vector po¬ 
tential can be disregarded i.e. A{1^) ~ Aq, one finds 
= EQ.l^{a + d’l') with ifo(^) = wpAg. This cor¬ 
responds to the electric-dipole approximation frequently 
used in atomic Cavity QED. In the opposite limit of a 
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small enough conductor or a perfect screening of Ai^) 
inside this conductor, one can use a constant in¬ 

side this conductor. This is reminiscent of the capacitive 
network models frequently used to describe supercon¬ 
ducting charge qubits in Circuit QED. Note that we do 
not assume that the ensemble T of the tunneling charges 
is neutral. 


C. Tunneling model and coarse graining of 
nanocircuits 

In this section, we consider a nanocircuit gathering 
quantum dots and reservoirs connected by tunnel junc¬ 
tions. The effect of time-dependent classical fields on 
tunneling through quantum dots has been studied since 
the 1990’s2^i^. Mesoscopic Circuit QED offers the op¬ 
portunity to study the interaction between microwave 
photons and quantum dot circuits from a different per¬ 
spective. Quantum dot circuits are often described in 
terms of a tunneling model in which the dots and reser¬ 
voirs contain quasi-localized orbitals coupled by tunnel 
elements^. This decomposition is instrumental in or¬ 
der to particularize the treatment of the reservoirs and 
account for the irreversibility of transport processes or 
other reservoir-induced damping effects. Importantly, 
the tunnel couplings arise from an overlap between neigh¬ 
boring orbitals. This means that these orbitals cannot be 
considered as orthogonal, except in the limit of weak tun¬ 
nel couplings (see for instance^^). However, in practice, it 
may be necessary to take into account a strong inter-dot 
tunneling. For example, in the case of a DQD, there can 
exist a strong splitting between bonding and antibond¬ 
ing states, which can be resonant with the cavityi^"— 
Depending on the absolute energy scale of the DQD con¬ 
finement potential, it might be necessary to go to the 
limit of a non-perturbative inter-dot tunnel coupling to 
obtain such a splitting. In order to circumvent the diffi¬ 
culty related to the non-orthogonality of the nanocircuit 
orbitals, we divide the nanocircuit into the ensemble T) 
of tunnel-coupled quantum dots from one side, and the 
different reservoirs with an index Rp with p £ {!,..., N}, 
from the other side. We use an index o G {V, Rn} 

to denote these different elements, and an index j to 
denote the different orbitals in a given element o. We 
assume that the low energy spectrum for the whole en¬ 
semble T) of the quantum dots can be determined in the 
absence of the reservoirs as the discrete spectrum of a 
potential profile with one or several potential wells. This 
procedure can be performed for an arbitrary inter-dot 
tunneling strength and leads to exactly orthogonal or¬ 
bitals, so that cx)j'} = 5jjr. One can also define 
exactly orthogonal orbitals in the isolated reservoir Rp, 
with p G {!,..., A^}, so that } = SpySjjr. 

In a second step, one can evaluate the tunnel coupling 
between the orbitals j and j' of V and Rp from 
the overlap between their wavefunctions. This gives 
= 0 only at lowest order in tunneling^. 


However, using a weak tunneling between 7) and Rp is not 
a severe restriction since the large density of states in the 
reservoirs can compensate for the smallness of the tun¬ 
nel coupling elements, and lead to a large tunneling rate 
compatible with the Hondo effect, for instance. To sum¬ 
marize, in the absence of the cavity, the tunnel Hamilto¬ 
nian writes 

He = ’^^SojchjCoj + ^ ^ + H.c.) (10) 

OJ o^o',j,j' 

where the creation operator corresponds to an or¬ 
bital j with energy Soj and wavefunction (poj{r) mainly 
localized inside element o G {V, Rn}- One can use 
{cljjCo'j'} = So,oAjj> at lowest order in the dot/reservoir 
tunnel couplings. 

We now reexpress the photonic pseudo-potential 
scheme of the previous section in the framework of the 
tunneling model. For this purpose, one needs to de¬ 
compose the field operator ^1 (r) associated to tunneling 
charges on the nanocircuit orbital states. At lowest order 
in the dot/reservoir tunneling, one can use^ 

o,j 

Hence, Hamiltonian ([H) directly gives 

=He + HhTui + knt (« + a+) + (12) 

with 

hint = '^^O'ojcljCoj + 'y ( -\- H.c.) . 

0,3 Oj^o'j' 

(13) 

Note that in principle, the term in from Eq.® renor¬ 
malizes Eoj, toj,o'j' and the Coulomb interaction term 
H^^i, but this is not essential for the physics we discuss 
below. 

Since the photonic pseudo-potential modifies the po¬ 
tential landscape seen by the tunneling charges, it can 
naturally affect all the parameters in the tunneling 
model. First, cavity photons shift the orbital energy Soj, 
with a coupling coefficient 

= -ejdr^ V±{f) . (14) 

In general, aoj strongly depends on the indexes o and 
j due to the space dependence of V_\_(r). This makes 
the cavity-induced orbital energy shifts particularly rel¬ 
evant experimentally. For a standard metallic reservoir, 
it is reasonable to disregard the dependence of aoj on 
the orbital index j, i.e. aoj — Oo, because the proper¬ 
ties of the electronic wavefunctions in the reservoir can 
be considered as constant near the Fermi energy^. In 
this framework, a behavior similar to the capacitive cou¬ 
pling of the cavity central conductor to the reservoirs 
is recovered. This type of reservoir/cavity coupling en¬ 
ables one to interpret data obtained for a quantum dot 



6 


coupled to two normal metal leads^d^. In contrast, aoj 
may strongly depend on j for o = V, as illustrated for 
instance by the example of a DQD with an asymmetric 

The Hamiltonian Herb also contains coupling terms in 

loo,o'r = -efdr^ip*j(^(po'f(r^V±(r^ (15) 

with oj ^ o'j'. These terms include photon-induced tun¬ 
neling terms in with o o', and photon-induced 

transitions internal to o, in In principle, all the 

above cavity/nanocircuit coupling terms can coexist. In 
practice, we expect the terms in aoj to be dominant 
over the tunnel terms in joj,o'j' with o o', which in¬ 
volve a weak overlap between wavefunctions. Neverthe¬ 
less, as shown in the next section, the signatures of the 
photon-induced tunneling terms can be boosted by the 
large number of states they affect in the reservoirs. If 
Vj_(r) varies slowly inside T), the photon-induced transi¬ 
tions internal to V are negligible, i.e. — 0, since 

l^PTij) and \ipvj') are orthogonal by definition. However, 
these transitions become possible if the cavity field is 
weakly screened ~ Eo.l^{d + al)), or if a pro¬ 

trusion of a cavity conductor is placed close to one of 
the quantum dots in order to reinforce the spatial vari¬ 
ations of Vj_ Transitions in 'yT>j,Vj' can be 

particularly interesting in case of real or artificial spin or¬ 
bit coupling which mixes spin and orbital states. In this 
case, the 'yx’j.'Dj' terms can correspond to spin-transitions 
inside a single quantum dot^ir— or a DQD^“— . In the 
reservoirs, the terms are probably always non¬ 

resonant due to strong energy relaxation. Hence, we will 
assume that these terms do not need to be treated explic¬ 
itly because they can be included in a renormalization of 
the states \ipR^jy 

In practice, the potential Vj_{f) can be evaluated nu¬ 
merically by removing all nano conductors from the de¬ 
vice and simulating the microwave electromagnetic field 
in the cavity in the presence of the metallic gates, sources 
and drains of the nanocircuit. This can be done using 
standard microwave simulation tools, which can take into 
account imperfect metals. It enables a realistic evalua¬ 
tion of the elements ttoj and Even if V±{r) is not 

calculated numerically but replaced by a phenomenolog¬ 
ical expression, the expressions m and (fT5)l remain in¬ 
teresting because they set constraints between the differ¬ 
ent aoj and 7oj,o'y, which all depend on the same Vj_(T) 
profile. 

We do not give details on the Coulomb interaction term 
^^oui which stems directly from Eq.Q. If Coulomb 
blockade in the nanocircuit is already strong in the ab¬ 
sence of the cavity, the effect of the cavity on interac¬ 
tions may be disregarded or treated perturbatively. In 
the next section, we discuss another situation where we 
assume that the effects of H^coui negligible due to the 
large tunnel rate between a quantum dot and a normal 
metal reservoir. In the general case. Coulomb interac¬ 
tions between the tunneling charges can lead to a large 
variety of effects, which we will not discuss in this work. 


To conclude, we obtain a decomposition of the system 
Hamiltonian in terms of nanocircuit elements connected 
by tunnel couplings. While this is reminiscent of the 
coarse graining description of superconducting microcir¬ 
cuits frequently used in Circuit QED, a full analogy is 
not possible due to the presence of the orbital degree 
of freedom in the Mesoscopic QED case. Due to this 
orbital degree of freedom, we obtain a large variety of 
electron/photon coupling terms. In this section, we have 
discussed the case of quantum dot circuits which raises 
most experimental efforts so far. However, the above 
approach can be used for any other type of system in 
which a decomposition in terms of a tunneling model is 
relevant. This includes quantum point contacts, molecu¬ 
lar circuits, metallic tunnel junctions, and hybrid super¬ 
conducting systems enclosing Andreev bound states or 
Majorana bound states, for instance. This last case is 
discussed in more details in section HT^ 


D. Discussion on the theoretical context 

Many theoretical works on Mesoscopic QED use as a 
starting basis a tunnel Hamiltonian where all parameters 
(orbital energies, tunnel rates,...) are perturbed by cor¬ 
rections proportional to the cavity electric field. These 
empirical approaches lack of a formal justification and 
can lead to unphysical predictions. Our work provides 
a rigorous framework for the use of tunnelling models 
with a linear light/matter interaction (the superconduct¬ 
ing case will be discussed explicitly in section HT^ . 

As already mentioned above, our approach proceeds 
along a spirit similar to the PZW transformation, but 
with modifications necessary to take into account the 
specificities of Mesoscopic QED, among which a strong 
space dependence of the cavity electric field on the scale 
of a nanocircuit, and boundary conditions provided by 
cavity conductors and voltage-biased DC electrostatic 
gates. Interestingly, the effect of boundary conditions 
provided by grounded conductors on the PZW transfor¬ 
mation has been discussed recently^. However, this Ref. 
considers only neutral atomic systems affected by locally 
constant cavity fields (dipolar coupling limit). 


E. Mesoscopic QED with superconducting 
nanocircuits 

1. Minimal coupling scheme with superconducting 
nanocircuits 

In the above sections, superconductivity was not ex¬ 
plicitly taken into account. It is important to discuss the 
generalization of our approach to the case of supercon¬ 
ducting hybrid nanocircuits, considering the present in¬ 
terest for Majorana fermions^^ or Andreev bound states 
for instance^. For simplicity, we consider below stan¬ 
dard BCS superconducting correlations characterized in 
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the mean field approach by a superconducting gap A(r^ 
coupling electrons and holes with opposite spins. The 
gauge invariant minimal coupling Hamiltonian of Eq. CD 
can be generalized as 

Hfot = Htot + f d^r + H.c.) 

(16) 

with '0|.(r) the field operator for tunneling electrons with 
spin cr € {tii} and 

$(r) = e(d — d^)V^{l^)/fujjQ . (17) 


with a spinorial wavefunction 

Wr.{z) =* [Ud;{z),ul{z),vl{z),-V^{z)] 

fulfilling 

heff(f)x„(f) = EnXn{r) (20) 

with Xn{,z) =* [u”(z), m"(z), —(z), —^ 1 ( 2 )] and 


heffif) 


hrir) A{r)ao 
A*(r)(To -c^yhr{r)(Jy _ 


( 21 ) 


The term Htot is the straightforward generalization of 
Eq.lJT]) to the spin dependent case, i.e. the single electron 
term can have a structure in spin-space, due to 

magnetic fields and/or spin-orbit coupling, for instance. 
Hence, V’(f^) must be replaced by *[^/-|-(r), i/)|(r)] in the 
first term of Eq. H]) and by'0-|-(r)-|-'(/’4,(r) in Eq.([3]). Due to 
gauge invariance, h'f{r) can involve the electron momen¬ 
tum through —iti^f+eA (r) only. Importantly, we define 
A(r) as a gauge invariant quantity, so that the phase pa¬ 
rameter <i> ensures the gauge invariance the pairing term 
of Htot- The form used above for $ is valid provided 
photon-induced magnetic effects can be disregarded, in 
agreement with the approach of section Hi HI 


Above, (To and ay are the identity and the second Pauli 
matrices in spin space. From the definition of c)^, one 
gets 


Htot = X! -I- huod^d 

n 

+ ('^nn'ClCn' + (d + O^) 

n,n' 

+ 12 + H-C.) (d + a^') 

n,n' 


with 


= -ejd^r w),(T^)14(T^)t^w„/(1^) , 


( 22 ) 


(23) 


2. Photonic pseudo-potential scheme with superconducting 
nanocircuits 

Using Eqs.® and CD with {r)'\j)(f) replaced by 
Scr'0j.(r)o.'0(r), Hamiltonian (fTHl) is transformed into 

H^ot = Htot + j d^r (|A(r)^|(70^j(f) -t H.c.) (18) 

where Htot is again the generalization of Eq.® to the 
spin dependent case. Importantly, Ea. dTBll shows that in 
the photonic pseudo-potential scheme, the phase of the 
gap term is not affected by photons anymore. This rep¬ 
resents one more advantage of this scheme. This result 
remains valid in the case of d-wave or p-wave supercon¬ 
ducting correlations^^. 


3. Mesoscopic QED with Bogoliubov-De Gennes equations 

In the case where electron-electron interactions can 
be disregarded, Bogoliubov-De Gennes equations are a 
widely used approach^, which enables a diagonaliza- 
tion of the circuit Hamiltonian in terms of quasiparticles 
with creation operator c)j. One can define quasiparticles 
modes in the absence of the cavity, by using 

ct = Jd^r (19) 


= -ejd^r Ul(I^) (|u]^(z)f + \v^z)\^) , (24) 

= -ejd^r -b iTj^X,(7^)/2 (25) 

and Tx, Ty , Tz Pauli matrices in Nambu space. 

4- Mesoscopic QED with Majorana fermions 

In this section, we discuss the case of a nanocircuit en¬ 
closing Majorana bound states, coupled to a cavity. In 
principle, these self-adjoint bound states can appear at 
the interface between the topological and non-topological 
superconducting phases of a nanostructure^izSS. Al¬ 
though an indirect Majorana bound state/cavity cou¬ 
pling mediated by a superconducting quantum bit or 
a Josephson junction has been considered in many 
works^^— , a direct coupling can also exist^. Using 
our formalism, we can provide a general Hamiltonian de¬ 
scribing such a situation. The direct Majorana bound 
state/cavity coupling occurs due an overlap between 
neighboring Majorana bound states, caused by the finite 
size of the nanocircuit. This overlap is naturally depen¬ 
dent on the photonic pseudo-potential. Proposals to ob¬ 
tain Majorana bound states in condensed matter systems 
rely on the use of superconducting elements, which calls 
for the use of the results of section llT^ For simplicity, we 
disregard Coulomb interactions and use the framework of 








the Bogoliubov-De Gennes equations of section lIlE3l To 
remain general, we do not specify the details of the sys¬ 
tem. We consider an ensemble of Majorana bound states 
with creation operators ^ = CM,n and wavefunction 
localized inside a nanoconductor. In general, 
the nanoconductor also encloses a continuum of ordinary 
fermionic states „ with wavefunction (ps,nC^) located 
above an energy gap. We assume that the nanoconduc¬ 
tor is tunnel coupled to a fermionic reservoir with orbital 
states Cjj „ with wavefunctions ifR^nC^)- Such a reser¬ 
voir can be used to measure the conductance through the 
Majorana systeni^. In the limit of weak tunneling, one 
can reexpress Eq. by using the ensemble of the oper¬ 
ators „ with o G {M,R,S}, following similar consid¬ 
erations as for the tunneling model of section III Cl This 
gives: 


ffS _ 


E 

n<.n' 


Cm ,n Cm,?! 


E 


^n^o,n^o,n 


+ CmA^u^u'Cr^^, - Cr^u') + 

n, n' 

+ E (^In.o'n'C'l.nC'o.n' +7l^^,o'n') 

o, o'g{M,S,R} 

n,n' 

+ E {lon,o'u'ClnCU, + H.C.) ^ + a) • 

o,o'g{M,S,R} 

n,n' 

( 26 ) 


Above, the coefficients 7 ^)) have expressions similar 

to the Vnn' coefficients of Eqs. (j231l25l) with Wn{n')(^) 
replaced by iPo{o'),n{n')(^)- We use energies for the 
orbitals in the reservoir (o = R) and the nanoconduc¬ 
tor continuum of states (o = S), and a tunnel coupling 
parameter On^n' between overlapping Majorana bound 
states. The Majorana bound states and the reservoir 
are connected by tunnel coupling terms in ■ For sim¬ 
plicity, we have omitted the DC tunnel coupling between 
the reservoir and the nanoconductor continuum of states. 
Photons can shift all the above mentioned parameters, 
due to the two last lines of Eo. d^ . In particular, pho¬ 
tons can modify the tunnel coupling between neighboring 
Majorana bound states, due to the terms in 7 ^^^ , as 

discussed in Ref^. Note that in Hamiltonian (I^Sl) . there 
is no DC coupling between the Majorana bound states 
and the nanoconductor continuum of states, which are 
assumed to be orthogonal by construction of CM,n and 
Cs,n'- However, there can be photon-induced transitions 
between the Majorana bound states and the nanocircuit 
continuum of states due to terms in 7 ^^^^ and 
Effects related to these last transitions have been dis¬ 
cussed in^. Our approach provides a general justifica¬ 
tion for such terms. 


III. EXAMPLE: QUANTUM RC CIRCUIT IN A 
CAVITY 

A. Hamiltonian of the system 

We now give an example of a specific prediction given 
by Hamiltonian m- We consider a cavity coupled to a 
” quantum RC circuit” made out of a single quantum dot 
tunnel contacted to a single normal metal reservoir. We 
assume that the dot has a large intrinsic level spacing 
so that a single dot orbital with energy Sd needs to be 
considered. The dot is capacitively coupled to a DC gate 
electrode, which enables one to tune Sd- We study the 
effects of the dot on the photonic mode with frequency 
ujQj2TT. Following the previous section, the total Hamil¬ 
tonian of the system is: 

HRC+cav = He -|- hinti& -f -f hujQa^a + Hbath + Hrr . 

(27) 

The Hamiltonian 

He =£dc\cd + jdeNQv{£)(ec?{e)c{£)+tc\c{e)+t*c'{e)cd) 

( 28 ) 

of the dot and its reservoir involves the lead density 
of states Nqv{£). For later use we have introduced a 
cutoff v{£) = /{D^ + £^). The terms in t describes 

tunneling between the dot and the reservoir, whose or¬ 
bitals correspond to fermionic creation operators and 
c’l’(e), respectively, with {cd, = 1 and {c(e), c^(e')} = 
(5(6: — £')/NoI'{£). One can use extra bosonic modes b^e) 
with [&(e), 61(6')] = S{e — e')/''^o to describe the intrinsic 
damping of the cavity, i.e. 

Hbath = f d£no{eb\e)b{e) + Tb\e)d + T*d^b{e)) . (29) 

The coupling between the quantum dot circuit and the 
cavity is described by the term hint{d -I- at) with 

hint = ac\cd -b Jd£Noiy(£)(jclc(£) + 7*ct(e)cd) . (30) 

It is not necessary to write explicitly the coupling be¬ 
tween the chemical potential of the reservoir and the cav¬ 
ity because since there is a single reservoir, it is only the 
difference between the potentials of the dot and the reser¬ 
voir which matters. Note that the terminology ’’quantum 
RC circuit” mainly refers to the fact that no DC current 
can flow through the device and dissipation is due to 
transport through a single tunnel barrier. However, in 
the general case, the concept of a local capacitor is not 
sufficient to describe the coupling of the electromagnetic 
field to the device, since the photonic pseudo potential 
modifies £d and t simultaneously. In practice, one can 
measure the response of the cavity to a microwave drive 
with frequency lorrAt:, which can be described by 

Hrf = £RFC‘^^^*a + . (31) 

Above, we have disregarded Coulomb interactions be¬ 
tween the tunneling charges (including the contribution 
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from the ftnoQ term arising from the unitary transfor¬ 
mation of section II.B), because we assume that these in¬ 
teractions are smaller than the tunnel rate V = tt |t|^ Nq 
between the dot and the reservoir. Such a simplifica¬ 
tion is frequently performed for interpreting standard 
current measurements in open quantum dot circuit a^^'^^ . 
In practice, one can check the relevance of this hypothe¬ 
sis by measuring for instance the AC current through the 
dot versus £d- 



FIG. 2: Scheme of the various contributions to Gf^. ^{t). 
The wavy lines correspond to photonic propagators, the sim¬ 
ple dashed lines to bare electronic propagators in the nor¬ 
mal metal reservoir, and the double full lines to electronic 
propagators in the quantum dot dressed by the dot/reservoir 
tunneling processes. 


with Go(a;) = (w — wq + *Ao)“^ the free cavity Green’s 
function. In the simplest model where r and rig are in¬ 
dependent of energy, the constant Aq = tt jrl^no/fi cor¬ 
responds to pure cavity damping. In the following, we 
consider the limit Aq ujq which corresponds to stan¬ 
dard cavities. In the limit ujrf — wq, Aq <C F, one can 
use a random phase approximation which leads to 


Ga,at(w_R_F) ~ (^UJRF — WQ + *Ao — H 


(33) 




the electronic re- 


spouse function calculated in the absence of the cavity. 
In the framework of a non-interacting diagrammatic cal¬ 
culation, ^(w) can be decomposed into the eleven 

contributions represented in Fig. 2. Our calculation is 
not perturbative in t since we use an exact expression 
for the dressed propagator of electrons on the dot, cor¬ 
responding to the double full lines in Fig.2. There are 
nine types of diagrams corresponding to the possible or¬ 
dered pairs formed with the interaction constants a, 7 , 
and 7 *. The diagrams of Fig.2a correspond to pure con¬ 
tributions from photon-induced orbital shifts in . The 
diagrams of Figs. 2b and c correspond to contributions 
from photon-induced tunneling in I 7 I , 7 ^ or 7 *^. Note 
that the contributions in I 7 I have been separated into 
the diagrams of Fig. 2b which diverge logarithmically 
for D large and the bottom diagrams in Fig.2c which are 
regular. There also exist interferences between photon- 
induced orbital shifts and photon-induced tunneling, as 
shown by the diagrams of Fig. 2d which depend on ay 
or ay*. 

The constants F and D set the scale of varia¬ 
tions of Gf I (uj). Therefore, in the limit wn <C 
T,D, one can replace “ 

Eg. (1551) . This gives a cavity frequency pull Awg = 
Ghi„A „X^ = 0 )/?i^ and a cavity damping pull AAg = 

In the limit F <C H 




-{uo/ih?) 

and at zero temperature, one obtains Awg = Awq 
where 


■ Awg 


B. Cavity frequency pull and damping pull 


Our purpose is to identify signatures of the terms in 
7 , which could become visible in the cavity response, 
considering the recent progress on the Mesoscopic QED 
technology22i2i. In the framework of the linear response 
theory, the cavity electromagnetic response to Hd is 
determined by the Green’s function^ Gd^i'^RF) with 

= 0)]) and Gj^j^{uj) = 
One can use the exact relation 

Ga,at(w) = Go(w) + ft-2Go(a;)G;,^„^;,^„^(a;)Go(a;) (32) 




(34) 


(Ftt^Aq Re[t7*]^F -|- 47ra Re[t7*]A’o£:d — a^F) 
stems from the diagrams of Fig.2,a, c and d, and 


F2 ^2 p 

D\og[ + \ed\ (2arctan[—] - tt) 


(35) 


stems from the diagrams of Fig.2,b. Note that the cut¬ 
off F^s) on the reservoir density of states is necessary to 
avoid a divergence of AuJq. The cavity damping pull is 
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given by 

^^0 = 2\2 “ 2TrRe[t-f*]Noedf ■ (36) 

+ ^d) 

Interestingly, Awq and AAg vanish for F, whereas 

Awq vanishes only for Sd^ D. Above, we have implicitly 
included the spin degree of freedom in the model. In the 
spinless case, Awq and AAq should be divided by a factor 

2 . 

From Eas. (l34ll3Sl) . at low temperatures, when the cav¬ 
ity has a frequency smaller than the dot/reservoir tunnel 
rate F and the reservoir bandwidth D, and when 7 = 0 , 
the cavity frequency pull and the damping pull fulfill the 
relation 0 = 7 r /2 with 

a^AAo , , 

ojo{hAuJo)^ 

This property is a manifestation of the Korringa-Shiba 
relational, whose universality is still actively discussed 
for quantum RC circuitslSr— . Remarkably, we find that 
this universality is broken for 7 finite. Hence, one can 
identify the existence of the 7 coupling by comparing the 
AuJo{ed) and VAAo(ed) curves (Fig. 3) or studying 0 
(Fig. 4). Here, we present results for 7 / 0 ; <C 1 because 
from Eqs. (nam), this is the most probable regime of 
parameters. When 7 = 0 , Awq and -x/AAo show iden¬ 
tical variations with e^, and are both even functions of 
Ed, while 0 = 7 r/ 2 . When 7 0, VAAo presents a res¬ 

onance with Ed which is wider than that of Awq (Eig.3, 
middle panel). A good resolution on the tails of the Alcq 
and a/AAq resonances could also reveal that Aujq and 
AAq are not even with Ed in the general case, except 
if arg[t 7 *] = 7 r /2 (Fig.3, bottom panel). In principle, 
even a relatively small 7 can be detected. Indeed, for 
7 / 0 ; ~ 10“'*, 0 already significantly deviates from 7 r /2 
at Ed 0 where the signals Awq and -^/AAq are maximal 
(Fig.4, pink curves). This is because the effect of the 
photon-induced tunneling term is boosted by the large 
number of states involved in the reservoirs. To conclude, 
in order to reveal the photon-induced tunneling terms, 
one promising possibility is to compare the cavity fre¬ 
quency pull and damping pull caused by a quantum dot 
with a single normal metal reservoir. Even a relatively 
weak photon-induced tunneling can affect these quanti¬ 
ties because the large number of reservoir states rein¬ 
forces the tunnel effect. 



-20 -10 0 10 20 

sjr 

a 


FIG. 3: Reduced Acjq and VAAq versus Ed for different values 
of 7 . We have used t/F = 0.001, D = 20r, and pulsation 
scales Ho = ‘lo? l-KVh and 17 1 = Anoflo/F. 



-20 -10 0 10 20 
syr 


FIG. 4: Ratio 0 versus Sd for different values of 7 . We have 
used the same values of t, D and f7o as in Fig.3 


C. Experimental state of the art and discussion 

In principle, circuit-QED should enable an accurate 
determination of 0 thanks to the use of resonant tech¬ 
niques. The fabrication of a quantum dot with a single 
normal metal reservoir inside a coplanar cavity is fully 
accessible with present techniques. However, so far, the 
closest case studied experimentally is a quantum dot with 
two normal metal reservoirs, which leads to more compli¬ 
cated physics due to an asymmetric coupling of the two 


reservoirs chemical potentials to the cavity2ii^. Different 
types of quantum RC circuits based on 2-DEGs but also 
carbon nanotubes or semiconducting nanowires deserve 
an investigation. The ratio 7 / 0 ; should be strongly de¬ 
pendent on the type of nanoconductor and the contacts 
configuration used. 

The physics discussed in this section is related to the 
problem of a quantum RC circuit which is not coupled to 
a cavity but directly excited by a classical microwave field 
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with amplitude urf and frequency a;/jF/27r, which corre¬ 
sponds to a term in hintaRF cos{ujRFt). At low frequen¬ 
cies and T = 0, it has been predicted^^— that a spinless 
quantum RC circuit excited classically should have the 
same admittance as an RC circuit, with a universal resis¬ 
tance Rrei = /i/2e^. However, this result was obtained in 
a purely capacitive coupling scheme i.e. 7 = 0. Assuming 
that the admittance of the circuit is determined from the 
current / = edfid/dt with hd = c^Cd, this value of Rrei 
can be recovered by calculating the charge susceptibility 
Gnd,nd(,^RF), he. the diagram of Fig. 2.a, which leads to 
Rrei = Qhje'^ = hl2e^ because 7 = 0. However, in the 
general case, one should calculate the current response of 
the quantum RC circuit from (^.^^{i^rf) because the 
oscillating drive can also modify the tunnel barrier trans¬ 
parency. Only the diagram of Fig. 2.a and half of the 
diagrams of Fig. 2d should contribute to Rrei, whereas 
the diagrams of Fig.2.b leading to the logarithmic term 
of Eg. ((551) should not contribute. Therefore, Rrei should 
show a behavior qualitatively different from 0. We will 
not discuss the detailed behavior of Rrei here. We will 
simply point out that for 7 0 , Rrei is not universal, 

and furthermore, Rrei and 0 are not trivially related. 
Interestingly, the universal value Rrei = hj^e^ expected 
for a spinless system with 7 = 0 has been checked experi¬ 
mentally in the limit T —>■ 0 for a 2-DEG based circuit^^. 
However, this measurement is not incompatible with a 
finite 7 , due to its uncertainty of the order of ± 20 % 2 ^. 

Note that since Rrei and 0 are qualitatively different 
signals, it can be interesting to measure both in the same 
experiment in order to obtain more information on the 
system. More generally, the study of tunneling effects can 
benefit from a simultaneous measurement of the cavity 
response and the current through the nanocircuit. Such 
a joint measurement has already been realized in the case 
of a DQD. The cavity dispersive shift is directly sensitive 
to the population imbalance between the DQD bound¬ 
ing and antibounding states whereas the current though 
the DQD corresponds to a more complex combination 
of state populations and tunnel rates. The simultaneous 
measurement of the two signals gives stronger constraints 
to determine the system parametersii. Another possibil¬ 
ity to gain more information on Mesoscopic QED systems 
might be to measure cross-correlations between the elec¬ 
tronic current and the cavity output field. 

Importantly, the example of the quantum RC circuit 
shows that dissipation of the cavity photons has to be 
considered carefully, since an open quantum dot circuit 
with a single contact can already induce photon dissipa¬ 
tion. On a more general level, photonic dissipation in 
Mesoscopic QED is a very rich problem. Another inter¬ 
esting possibility would be to prepare non-trivial quan¬ 
tum photonic states by combining the effects of photonic 
and electronic dissipation with the non-linearity of quan¬ 
tum dot circuits, by analogy with the methods devel- 
opped for two level systemsS^i^l. This so-called ’’reser¬ 
voir engineering” could benefit from the specificities of 
nanocircuits like for instance the use of reservoirs with 


specific electronic orders or finite bias electronic trans¬ 
port . 


IV. CONCLUSION 

Mesoscopic QED has common ingredients with atomic 
Cavity QED, like the relevance of an orbital degree of 
freedom, which does not exist for superconducting qubits. 
It has also common ingredients with standard Circuit 
QED, like the tunneling physics and the strong inhomo¬ 
geneities of the photonic modes. It is therefore necessary 
to develop a specific description which combines all these 
ingredients. In the case where photon induced magnetic 
effects can be disregarded, one can express the coupling 
between tunneling quasiparticles and photons in terms 
of a scalar photonic pseudo-potential. In the framework 
of a tunneling model, this leads to photon-induced or¬ 
bital energy shifts, which coexist with photon-induced 
tunneling terms and photon-induced local orbital transi¬ 
tions. To illustrate the richness of our approach, we have 
discussed the example of a cavity coupled to a quan¬ 
tum RC circuit, i.e. a single quantum dot coupled to a 
single normal metal reservoir. The photon-induced tun¬ 
neling terms between the dot and the reservoir induce a 
non-universal relation between the cavity frequency pull 
and the cavity damping pull, contrarily to what is ex¬ 
pected with purely capacitive coupling schemes at low 
temperatures^^— . This case represents only one exam¬ 
ple of use for our formalism. We have given explicit 
Hamiltonians for the cases of a multi-quantum dot cir¬ 
cuit in a cavity, and a superconducting nanostructure en¬ 
closing Majorana bound states coupled to a cavity. One 
could also consider cavities coupled to quantum point 
contacts, molecular circuits, and metallic tunnel junc¬ 
tions, for instance. 
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Saptsov, M. Trif, and J. J. Viennot. This work was 
financed by the ERC Starting grant CirQys, the EU 
EP7 project SE2ND[271554] and the ANR-NanoQuartet 
[ANR12BS1000701] (France). 


V. APPENDIX A: EFFECTIVE MODEL FOR 
MESOSCOPIC QED DEVICES 

A. Definition of the model 

In the main text, Eqs O-® have been introduced 
on the basis of a physical discussion and gauge invari¬ 
ance considerations. It is also instructive to use an effec¬ 
tive model which separates spatially the tunneling and 
plasmonic parts of the nanocircuit. Note that this effec¬ 
tive model does not aim at calculating quantitatively the 
nanocircuit/cavity couplings. It rather aims at providing 
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FIG. 5: (a) Effective model for the circuit of Fig.l, where 
the nanocircuit reservoirs are decomposed into effective or¬ 
bital reservoirs (dark grey elements) and effective plasmonic 
reservoirs (nearby purple elements) (b) Ensemble C of the 
perfect conductors considered for the generation of the pho¬ 
tonic modes. The photonic field in yellow is inhomogeneous, 
as represented by darker yellow areas near the cavity protu- 
sions, and by the white (screening) areas between the small 
purple conductors (c) Effective orbital nanocircuit O. 


a full mathematical justification for the form of Eqs. O- 
(|4|) and a deeper insight on the underlying physics. For 
simplicity, we assume that the nanoconductors host only 
individual electronic orbital modes since their electronic 
density is generally very low. We assume that the plas¬ 
monic modes in the fermionic reservoirs (blue elements 
in Fig. 1) are ultra fast due to the absence of dynamical 
Coulomb blockade^. In this case, one can decompose 
each fermionic reservoir into a purely orbital conductor 
(dark grey elements in Fig.5.a) which contains no plas- 
mons and is tunnel coupled to the nanocircuit, and per¬ 
fect conductors (small purple elements in Fig.5.a) which 
can host plasmonic modes, but are only capacitively cou¬ 
pled to the nanoconductors. These perfect conductors, or 
’’effective plasmonic reservoirs”, can be placed such that 
the partial or full decrease of the photonic fields around 
the nanocircuit orbital states is reproduced at least qual¬ 
itatively. We denote with C the ensemble of the perfect 
conductors, which includes the cavity conductors, gate 
electrodes, and the effective plasmonic reservoirs (see Fig. 
5.b). Some of these conductors are voltage biased and 
some others are left floating, like for instance the central 
conductor in a coplanar cavity. The rest of the circuit 
represents an effective orbital nanocircuit O where tun¬ 
neling physics prevails (see Fig. 5.c). Note that in prin¬ 
ciple, to grant current conservation, the ’’orbital” reser¬ 
voirs have to be connected to voltage sources through 


wirings which host plasmonic modes, but one can as¬ 
sume that they are too far from the nanoconductors to 
have any significant influence on the nanocircuit/cavity 
coupling. 

For later use, we attribute to each conductor of the en¬ 
semble C an index i and we note Si and Vi the correspond¬ 
ing surface and volume. In the following treatment we 
make a distinction between the conductors i G F which 
are left floating with a constant total charge Qi, and the 
conductors I G B which are voltage-biased with a con¬ 
stant voltage Vi or grounded (V) = 0). The generators 
imposing the voltages V) can be omitted from the de¬ 
scription and included in electromagnetic boundary con¬ 
ditions since they are far from the cavity and nanocir¬ 
cuits. The volume Vo complementary to ^iVi hosts a 
charge distribution t) = 5{^ —q a) of par¬ 

ticles a with charges Ca at positions if a, which belong 
to the nanocircuit O. We also define the corresponding 
current distribution j = X]aGC>®“ 9 ~ ^a)- 

B. Decomposition of fields in harmonic, 

longitudinal, and transverse components. 

In Cavity QFD, the treatment of the boundary condi¬ 
tions provided by the cavity conductors is usually omit¬ 
ted, because the atoms are very far from the cavity mir¬ 
rors. Therefore, a procedure based on a spatial Fourier 
transformation can be used to separate the transverse 
electric field which has no divergence, from the 

longitudinal electric field i^||, which has no rotational^. 
It is found that the value of if || is instantaneously im¬ 
posed by the charge distribution p{l^, t) associated to the 
atomic system. In contrast, the transverse fields and 
if = if ± correspond to propagating modes. In Meso¬ 
scopic QFD, it is necessary generalize this procedure to 
account for the presence of the perfect conductors of the 
ensemble C, including voltage-biased electrostatic gates. 
As already mentioned in the previous section, we use a 
charge distribution p{l^,t) which includes the charges 
from the O subsystem only. We will take into account 
the screening charges on the surface of the conductors of 
C through electromagnetic boundary conditions. 

Inside Vq, p{l^,t), j ("r^,!), and the total electric 
and magnetic fields and lf{l^,t) are related by 

Maxwell’s equations. To understand how these relations 
are affected by the screening charges lying on C, we in¬ 
troduce the Hodge decomposition of a generic field / as 
f = f\\ + fj- + harm, where /|| has a finite divergence 
but no rotational, f± has a finite rotational but no di¬ 
vergence, and fharm has non o^^FO . A similar decompo¬ 
sition was invoked very recently in the context of Cavity 
QFD^. However, in this Ref., the effects of the DC volt¬ 
age biases Vi or the floating charges Qi are not considered 
(i.e. T( = 0 is used for any i). Besides, the volume Vo 
outside the cavity conductors is assumed to be bounded. 
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This hypothesis is not possible in our case due to open¬ 
ings in the cavity planes, which are necessary to connect 
the nanocircuits to their outside voltage bias circuitry. 
Hence, we need to perform the Hodge decomposition in 
the fully general case. Below, we explain how to perform 
this task. 

It is convenient to introduce a potential ^harmC^) 
and a Green’s function G(l^, 1^') which fulfill the static 
Laplace equations: 


can check that the longitudinal components and j || 
are set instantaneously by since 

= p(7^,t)/eo (49) 

for G Vo- Furthermore, ^_l, and j ± follow an 
independent system of propagation equations, i.e. 

= (50) 


and 

^y^^harm(~^) — 0 

(38) 



(39) 

for (7^,1^') G Vo, with boundary conditions 



fg_d^r7^^harmC^)-n = -Qi for z e F , 

(40) 

and 

^harmC^) = for 1^ G Si and z S B 

(41) 


A = po7±(y,t) + (51) 

at 

for G Vo- The above equations have to be supple¬ 
mented with the boundary conditions 

/g.d^ri?||(_L)(l^,t).n = 0 for i S F , (52) 

.^.i?ll(j_)(T^,t) = 0 for G Si with (z, j) G (53) 
with £ij any trajectory connecting conductors i and j. 

C. Mesoscopic QED classical Lagrangian. 


d^r^^G(7^, 7^').n = 0 for z G F , 

(42) 

G(l^, 1^') = 0 for 7^ G Si and z G B . 

(43) 

Above, ft is the outward pointing unit vector perpendic¬ 
ular to Si- One can check that each of these two sets of 
Eqs. has a unique solution provided B is not empty^. 
From Eos. (I381I43D and the Maxwell Eos., one can check 
that: 


(44) 

4(^,0 = -/dV^^G(T^,T^')p(^',i) 

(45) 

and 


E±{'r^,t) = EC^A) - “ Eharmiv^) 

(46) 

while is still valid. We also find 

j = j \\ + j ± with 

7\\C^,t) = -('94||(l^,t)/at)/cVo 

(47) 

and 


h 

II 

1 

^■'1' 

(48) 


The harmonic component Eharm represents a static con¬ 
tribution to the electric field, set by the charges Qi or 
potentials Vi imposed on the conductors i- Thanks to 
the use of 0(1^, ^'), the definitions (I45M8I) take into ac¬ 
count the effect of the screening charges lying on C. Using 
the Maxwell equations and the expressions (I44II48I) . one 


At this stage, it is convenient to introduce scalar and 
vector potentials U and A in the Coulomb gauge i.e. 
7^.A{^,t) = 0- We also define the decomposition 
U{~^A) = + 4*/iarm(T^), such that 



(54) 

= -dAil^,t)/dt , 

(55) 

B — 7^ A A{^, t) . 

(56) 


A comparison between Eqs. (I39II43I) . and (I49I52I53I) leads 
to the identification 


C/||(T^,t) = (57) 

up to a global constant which one can disregard. Hence, 
A|| but also C/|j are instantaneously set by the charge dis¬ 
tribution p{l^,t)- In contrast, the transverse fields 
and B are determined by A which follows the propaga¬ 
tion equation 

□v(A(7^, t)) = -Aio7-L(^, t) (58) 

given by Ea. dCT) . with — (9^/c^9t^), and 

boundary conditions similar to Eas. (|5^ - (IM1) . In this 
framework, a proper classical Lagrangian is 


L = ^J2arncql - t) 

harmi~^ a) T ^ a ■ 4^ ( ^ a 


( 59 ) 


eo 


L (fr 

2 *^Vc) 


dt 


-c^ 
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with a € O in the sums. The Lagrangian equations de¬ 
riving from L correspond to Eq. ((551) and to the standard 
Newton equations for the motion of particles a. 


one can make a numerical microwave simulation using 
the real device geometry, with the nanoconductors omit¬ 
ted (see main text). 


D. Mesoscopic QED Hamiltonian 


VI. APPENDIX B: ADVANTAGES OF THE 
PHOTONIC PSEUDO-POTENTIAL SCHEME 


From Eq. (l59l) , one can check that Pa = ma q a + 
a,t) is the conjugate of a and n_L(T^,t) = 
—eoE±(l^,t) is the conjugate of Therefore 

one can express the Hamiltonian of the system as 

H L. This 

Hamiltonian can be quantized straightforwardly by re- 

placi^ "p 
and A. This gives 

I 


^ I? and ^ by operators pa, q a, H 


H = 


2m, 




(60) 


+ 9 «) + J2a^a^h 


2 J Vc) ^ 


(1 


^ I 

— 

H 

-f — 

\£o 


Mo 


arm\ Q ck) 


A(T^) 


In principle, the ensemble of the charges a G O includes 
tunneling electrons {a G T) but also ions and valence 
electrons from an underlying crystalline structure. In 
the main text, the two latter are assumed to be decou¬ 
pled from the cavity and are thus treated in a mean field 
approach. Note that the harmonic potential term in the 
above Hamiltonian is omitted by Ref.—. In our case, this 
term is crucial to account for the effect of DC electrostatic 
gates. 

Assuming that th^coupli^ between C and O is per¬ 
turbative, the fields A and H can be expressed in terms 
of photonic modes calculated in the absence of O, using 
for instance Ea. ((58)l with j ± = 0 We treat explicitly 
only one of these modes in Eq.(|5]). Screening charges 

on the nanocircuits reservoirs can strongly modify A. 
This is why, in the general case, it is not adequate to 
make a perturbative treatment between the photonic field 
generated by the empty cavity and the whole nanocir¬ 
cuit. Instead, the present work considers a perturba¬ 
tive treatment between the photonic field generated by 
the plasmonic modes in the whole Mesoscopic QED de¬ 
vice, and the effective orbital nanocircuit O. Similarly, in 
superconducting Circuit QED, the cavity modes can be 
strongly renormalized by the presence of a superconduct¬ 
ing quantum bit, and perturbation schemes must there¬ 
fore be defined carefullySS. 

In this Appendix, the renormalization of the bare cav¬ 
ity modes by the nanocircuit plasmons is described for¬ 
mally by supplementary boundary conditions imposed 
by effective plasmonic reservoirs. This enables a simple 
mathematical justification for the form of Eqs (II])-(l4]). In 

order to get realistic estimates of the fields A and H, 


The advantages of a formalism without the A^ terms 
have already been discussed thoroughly for Cavity QEDi. 
However, since the physical ingredients and measured 
quantities are different in Cavity and Mesoscopic QED, 
it is useful to discuss more specifically the advantages of 
the photonic pseudo-potential scheme of section HIE 21 
The first quantity which is generally measured in Meso¬ 
scopic QED experiments is the cavity frequency pull for 
a given cavity mode, which is modified by the nanocir¬ 
cuits coupled to the cavity. Since the minimal coupling 
scheme and the photonic pseudo-potential scheme are re¬ 
lated by a unitary transformation, they must of course 
predict the same cavity pull. It is nevertheless instructive 
to check this equivalence for a simple example. Here, we 
use a simpler approach than in section IV, because we 
do not need to particularize the effect of the nanocircuit 
reservoirs and we will not discuss the cavity damping 
pull. We consider a cavity coupled to a nanocircuit with 
unperturbed many electron states \ipj) which satisfy the 
general equation Hq \ipj) = Ej \ipj) with 


/ ^2 \ 2 

Hq = ^ -e$tot( q a)j + q a, q a') 

^ ( 61 ) 

and 4>tot( q a) = ^harmi q a) + Vconfi q a) ■ We discuss 
the limit where the nanocircuit and the cavity are not 
resonant i.e. Ej — Eji ^ hujQ for any j and j'. At second 
order in the cavity/nanocircuit coupling (i.e. V± or A), 
an elementary perturbation theory gives the cavity pull 


AHuJoj = 


2{nu:^f{E,-Ey)M^P 


(62) 


with 

I'fj') (‘ff I y±i^') I'fj) ■ (63) 

Note that Ahujoj depends on the state \ipj) occupied by 
the nanocircuit. This can be used, in principle, for a 
noninvasive readout of the nanocircuit state in the non¬ 
resonant regime^. Equation (I62|) can be obtained by 
using Hamiltonian ([5]) or equivalently Hamiltonian (fT|). 
In order to check the equivalence between the two ap¬ 
proaches, it is necessary to invoke the completeness re¬ 
lation \ipj) {ipj \ = I for the nanocircuit states. Inter¬ 
estingly, one can check that the A^ terms in Eq. © give 
a contribution 

E 2(E,-E,QA4““ 

j' ,a,Q!' 


(64) 
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to Afiwoj, which is not negligible for \Ej — Eji \ ^ Swq- 
We conclude that in the framework of Mesoscopic QED 
where the cavity frequency pull is a central quantity, it is 
particularly important to eliminate in a rigorous way the 
terms if one wants to use a linear light/matter cou¬ 
pling. The photonic pseudo-potential scheme completes 
such a task. 

The cavity frequency pull is not the only means to 
characterize the interactions between a nanocircuit and 
a cavity. A more general workout of Mesoscopic QED 
devices requires the knowledge of the matrix elements 
between the different nanocircuit eigenstates, generated 
by the cavity photons. Another argument in strong favor 
of the use of the photonic pseudo-potential scheme is the 
general relation 




ie 

2ma{Eji — Ej) 


{ipj\ p a-A{qa) + A{qa). p al^Pj’) ■ 

(65) 


It indicates that the single-photon coupling elements de¬ 
cay more quickly with [Eji — Ej) in the photonic pseudo¬ 
potential scheme. In practice one often uses a truncated 
space for the nanocircuit electronic states which is dif¬ 
ficult to handle or determine globally. This is what we 
do for instance in section IV, where we consider a sin¬ 
gle orbital quantum dot. Equation (1651) shows that in 
this case it is more accurate to use Eq.([5]) to predict the 
cavity behavior. 


Einally, the photonic pseudo-potential scheme involves 
only photon-independent superconducting gap terms (see 
section III El) . This can be another significant advantage 
considering the already rich structure of the nanocir¬ 
cuit/light coupling in the absence of superconductivity. 
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